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ABSTRACT 

The upper envelope of the amphtude of the VLBI visibihty function usually represents 
the most compact structural pattern of extragalactic radio sources, in particular, the 
"core-jet" morphologies. By fitting the envelope to a circular Gaussian model in ~3000 
parsec-scale core-jet structures, we find that the apparent angular size shows significant 
power-law dependence on the observing frequency (power index n = — 0.95±0.37). The 
dependence is likely to result from synchrotron self-absorption in the inhomogeneous 
jet and not the free-free absorption [n = —2.5), nor the simple scatter broadening 
(n < -2). 
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1 INTRODUCTION 



Radio-loud active galactic nuclei (AGNs) are a class of high- 
luminosity objects in the Universe. Their pc-scale jets can 
be well imaged by VLBI observations. The "angular size - 
frequency" ("0 - f") relation can be used for calibrate the 
apparent "angular size - redshift" ( "S - 2" ) relation in which 
the radio jets are taken as a "standard rod" applicable for 
testing different cosmological models (Gurvits 2003). The 
dependence of the size on frequency can provide impor- 
tant constraints for studying the internal physical properties 
(Lobanov 1998) and the evolution of jets. 



An extragalactic radio source generally has a power-law 
total flux density spectrum resulting from the synchrotron 
radiation. However, when imaged with VLBI at a milliarc- 
second resolution, the jet emission often exhibits more com- 
plicated spectral shapes, affected by opacity in the emitting 
material itself and in the surrounding medium. The interven- 
ing media along the line of sight scatter/absorb the emission 
and distort the image, in particular at low (< 1 GHz) fre- 
quencies. Thus, the apparent angular size usually depends 
on the observing frequency. The dependence can be written 
phenomenologically as a power-law function 8 oc in most 
cases. For instance, n — —2 for the compact radio source 
at the center of our Galaxy (Sgr A*) as a result of scatter 
broadening (Shen et al. 2005). 



Table 1. Summary of the used data set. 
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2 THE CHARACTERISTIC ANGULAR SIZE 

The characteristic angular size is a parameter describing the 
width of the intensity distribution of jet emission. Gener- 
ally, it can be defined either in the image plane or in the u-v 
(Fourier) plane. For instance, a measure of the characteristic 
size can be provided by the distance between the brightest 
component and the most distant one with the brightness at 
least 2 % of the peak (Kellermann 1993) in the image plane. 
It is also feasible to obtain the estimated angular size by fit- 
ting a circular Gaussian model to the visibility data (Gurvits 
1994) in the u-v plane. Although the two ways are equiva- 
lent in principle, the definitions in the u-v plane are more 
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Figure 1. An example for the definition of tiie cliaracteristic 
angular size. The top panel shows the original visibility data. The 
bottom panel shows the extracted upper envelope points and the 
fitting Gaussian curve. 



direct than those in the image plane, and free from sam- 
pUng effect. If the visibihty data have enough sensitivity, the 
"super-resolution" can be obtained for the barely resolved 
sources using a certain model (Lobanov 2005). Here we de- 
fine the characteristic angular size 9 by fitting the upper 
envelope of the visibility amplitudes to a circular Gaussian 
model (Pearson 1999): 



r(p) = exp 



41n2 



(1) 



where 6 is the full width to half-maximum intensity 
(FWHM) in the image plane in radian, p is u-v radius in 
wavelength and r(p) is the normalized amplitude. The char- 
acteristic angular size represents the most compact and the 
brightest part of the source. Fig. 1 gives an example for the 
definition of the characteristic angular size. 



3 DATA COLLECTION 

Table 1 gives a brief summary of our data set of 4000 objects 
compiled from several large VLBI surveys (see references 
in Prey 2006). Some sources have multi-frequency and/or 
multi-epoch observations. All the visibility data have been 
calibrated, self-calibrated, and saved as standard u-v PITS 
files. We developed several programs to automatically calcu- 
late the characteristic angular size according to our defini- 
tion using the IDL Astronomy Users Library. We extracted 
the upper envelope points from the original visibility data, 
fitted them by a circular Gaussian model, and finally fitted 
the "0 - v" data to a power-law function to determine the 
power index n for the multi-frequency observations. 



4 RESULTS 

Pig. 2 shows the statistical distributions of the characteristic 
angular sizes. The histograms show a similar shape peaking 
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Figure 2. The statistical distributions of the characteristic angu- 
lar sizes. Different colors represent different bands. Note that the 
angular size is expressed in units of the resolution of the VLBA 
array at each band. 



at the position of ~0.4i?, where 7? is the resolution of the 
VLBA array at each band and listed in Column 3, Table 1. 
The Ku-band distribution is consistent with the previous re- 
sults in the MO J AVE sample by Kovalev et al. (2005). They 
fitted the core by an elliptical Gaussian model in Difmap. 
Fig. 3 shows the distribution of the power index n, which has 
a mean value of —0.95 ± 0.37 for 2827 sources with at least 
two-frequency observations and —0.96 ±0.21 for 315 sources 
with at least four-frequency observations. The dependence 
agrees well with the results from fitting the core by an el- 
liptical Gaussian model in 167 objects (Jiang et al. 2001) . 
Note that we assume that the power index n is independent 
of redshift within the studied range of redshift (0 < z < 5). 



5 DISCUSSION 

The upper envelope statistically gives the size of the bright- 
est and most compact component. Among all the compo- 
nents in a pc-scale jet, the core/jet base usually has a fiat- 
ter spectrum and dominates the total flux density. Thus, the 
upper envelope gives the size of the radio core in most cases. 

The radio core generally has a flat spectrum and is 
located in the region where the optical depth r = 1. 
Konigl (1981) demonstrated that, due to the synchrotron 
self-absorption, the size of the core rcoro varies with the fre- 
quency V as: Tcore oc whcrc fcr depends on the shape 
of the electron energy spectrum and on the distribution of 
the magnetic fleld and particle density in an inhomogeneous 
jet, fcr = 1 in the case of the equipartition between jet parti- 
cle and magnetic field energy densities (Lobanov 1998). Our 
results agree well with this prediction. 

Free-free absorption by the plasma covering the jet also 
results in an optically thick core. However, free-free absorp- 
tion is only found in a few sources. Lobanov (1998) suggested 
that the shift of core position is rcorc oc //^^'^ for a spher- 
ical distribution of free- free absorption plasma Ac oc r~^, 
where r is the distance from the center and Ae is the elec- 
tron number density. So, the free-free absorption is not the 
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Figure 3. The statistical distribution of the power index n. The 
solid line represents the result from all the sources with at least 
two-frequency observations. The dotted line represents the result 
from all the sources with at least four-frequency observations. 



right mechanism to explain the statistical relation presented 
here. It might however happen in some sources. As far as 
the simple difFractive scattering models, e.g. the Gaussian 
screen model or the power-law model with the Kolmogorov 
spectrum, are concerned, they result in a dependence with 
n < — 2 (Thompson et al. 1986), which is fairly far from the 
peak (n ^ —0.9) of the obtained distribution. Thus, simple 
scatter broadening cannot explain the dependence. 
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6 CONCLUSIONS 

Based on the above discussion, we can draw the following 
conclusions. 

(i) The upper envelope of visibility amplitudes is a sta- 
tistically meaningful characteristic angular size for the radio 

"core" . 

(ii) The apparent "angular size - frequency" dependence 
can be statistically described by 6» oc ;y-o-9B±o.37 -^^ ^3000 
pc-scale jets. 

(iii) The observed frequency dependence can be explained 
as a result of synchrotron self-absorption in the iiihomoge- 
neous jet. The simple scatter broadening and free-free ab- 
sorption can be ruled out as dominant mechanisms respon- 
sible for the dependence. 



7 ACKNOWLEDGMENT 

J. Yang is grateful to the KNAW-CAS grant 07DP010. 
S. Frey acknowledges the OTKA K72515 grant. We made 
use of the collection of calibrated VLBI u-v data main- 
tained by L. Petrov at NASA Goddard Space Flight Cen- 
ter (GSFC) . This research has made use of NASA's Astro- 
physics Data System, NASA/IPAC Extragalactic Database 
(NED) and the United States Naval Observatory (USNO) 
Radio Reference Frame Image Database (RRFID). We 
thank A.L. Fey for providing the RRFID data. The National 



